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INTRODUCTION 

The  fluidization  technique  has  been  employed  widely  in  many  fields  of  the 
process  industry.  The  term  fluidization  generally  implies  that  solid  particles 
are  suspended  in  a  fluid  flowing  through  a  vessel.  The  mixture  of  particles 
and  fluid,  therefore,  forms  the  over-all  bed. 

The  advantages  and  disadvantages  of  the  fluidization  technique  in  com- 
parison with  the  conventional  packed  bed  reactors  can  be  found  elsewhere. 
(11,  28,  2$) ,  May  (15>)  in  his  recent  paper  pointed  out,  "Chemical  reactions 
with  fluidized  beds  have  shown  efficiencies  varying  all  the  way  from  'better 
than  fixed  beds'  to  as  low  as  $%   of  fixed  beds.  The  low  figure  shows  that  the 
bypassing  must  be  important ," 

Milhelm  and  Kwauk  (26)  first  observed  and  described  the  distinction 
between  particulate  fluidized  bed  and  aggregative  fluidized  bed.  In  the 
particulate  fluidization,  particles  are  discretely  separated  from  each  other 
and  a  mean  free  path  for  the  particles  seems  to  exist j  whereas,  the  particles 
in  the  aggregative  fluidization  are  present  in  the  bed  not  as  individual  units, 
but  as  local  aggregates.  The  particulate  fluidized  bed  resembles  closely  an 
ideal  expansion  of  a  fixed  bed,  whereas  an  aggregative  fluidized  bed  is  a 
heterogeneous  body  composed  of  two  separate  coexisting  phases  —  dense 
phase  and  bubble  phase  (dilute  phase) <>  A  close  analogy  can  be  found  between 
the  characteristics  of  a  fixed  bed  and  a  particulate  fluidized  bed,  but  it  is 
not  so  between  a  fixed  bed  and  an  aggregative  fluidized  bed.  The  differences 
are  chiefly  due  to  the  fact  thats 

a)  Particle  distribution  in  the  aggregative  bed  is  extremely  heterogeneous 

b)  There  is  no  definite  way  to  measure  the  height  of  bed  expansion. 


Literature  Survey 


During  the  past  years,   the  characteristics  of  aggregative  fluidized  beds 
have  been  studied  by  several  different  experimental  methods.     The  experimental 
methods  used  by  previous  workers  are  summarised  as  follows? 

1.  High  speed  photographic  studies  were  used  by  Mathcson  et.al.(llj), 
Furukawa  and  Gmae  (7)r  and  by  Leopoldo  Massimilla  and  J.  W.  Westwater   (13) • 
They  took  high  speed  photographs  of  bubbles  through  transparent  walls  of 
fluidized  columns  for  gas-solid  systems. 

2.  Shuster  and  Kisliak  (2U)   studied  the  uniformity  of  aggregative 
fluidized  beds  by  using  diaphram  plates  to  measure  pressure  gradients. 

3.  By  measuring  the  variation  of  the  electric  capacitance  between 
condenser  probes,  Morse  and  Ballou  (18),   and  J.  M.  Dotson  (5)  examined  the 
uniformity  of  aggregative  fluidization.     Later  Bakker  and  Heertjes  (1,   2) 
applied  a  similar  technique  to  measure  the  point  porosity  of  fluidized  beds 
at  various  locations. 

lu     Yasui  and  C.  N.  Johanson  (27)  devised  a  small  light  probe  to  study 
the  frequency,  thickness  and  the  rising  speed  of  bubbles  in  aggregative 
fluidized  beds. 

5.     The  X-ray  absorption  technique  made  it  possible  for  Grohse  (9)  to 
measure  the  over-all  density  of  the  aggregative  fluidized  bed  without  dis- 
turbing the  bed.     Baumgarten  and  Pigford  (3)   employed  a     #"-ray  absorption 
technique  in  the  study  of  bubble  size  and  bubble  frequency  in  aggregative 
fluidized  beds. 


*  A  rather  complete  literature  survey  concerning  fluidization  can  be  found 
in  the  Ph.D.  thesis  by  P.  T.  Shannon  (22). 


Some  of  these  investigators  have  proposed,  certain  mathematical  or 
empirical  model?  to  fit  their   experimental  data0     None  of  these  models  has 
yet  been  verified  as  acceptable  and  applicable  to  aggregative  fluidized  beds 
in  general-. 

Purpose 

The  major  purpose  of  the  present  work  was  to  investigate  the  axial  dis- 
tribution of  the  fluidized  bulk  density  in  an.  aggregative  fluidized  bed  by 
2T-ray  technique.     This  is  to  determine  qualitatively  and  quantitatively  the 
axial  density  profile  as  the  function  of  several  operational  variableSo 

EXPERIMENTAL 


Application  of    ?T-ray  Absorption  Technique 
to  a  Fluidized  Bed 


The  radiation  attenuation  method  of  measuring  density  of  a  twc-phase 
fluid  has  become  prevalent  during  the  past  years  (19,   3)»     This  method  of 
measuring  density  is  based  on  the  absorption  of  high  energy  radiation  from 
a  radioactive  source  which  can  be  detected  by  a  scintillation  counter  probec 
Two  of  the  characteristics  of     tf"~rays  are  their  high  energy  and  their  short 
wave  length.     When  it   is  applied  to  the  measurement   of  density,    it   is  essen- 
tially a  process  of  counting  the  number  and  indicating  the  kind  of  atoms  by 
which  the  radiation  is  being  absorbed.     The  attenuation  of     sT-rays  is  vir- 
tually independent  of  the  chemical  or  physical  state  of  the  atoms „ 


If  the  absorption  of  the  primary  radiation  by  air  (or  by  other  fluidizing 
gas)  is  negligible,  the  attenuation  of    ^f* -rays  shown  on  the  scintillation 
probe  is  due  to  the  absorption  of     a  -rays  by  the  thickness  of  the  walls  of 
the  fluidizing  column,,  the  constant  absorbing  media,  such  as  the  protective 
enclosure  of  the  scintillation  crystals  (Fig.   1),  and  the  particles  packed  or 
fluidizing  in  the  path  between  the  radioactive  source  and  the  scintillation 
counter  probe.     The    ^-rays  describe  an  essentially  parallel  path  in  passing 
from  the  source  to  the  scintillation  counter  probe.     The   #"~-rays  absorption 
by  a  homogeneous  material  of  density,     P  ,   and  thickness,  D,  can  be  described 
by  the  Lamberts -Beer's  Laws 

IsIne  ' 


-o 

or 


la(  I/I0)   -     -JJ^f  I  (1) 

The  absorption  coefficient  " U    »  varies  in  accordance  with  the  following 
relationships 

Mn  *  (VM)  tf~ 

where, 

Wa  is  avagadro»s  number 

M  is  the  atomic  weight,  and 

CT  Is  the  microscopic  absorption  cross  section. 

Therefore,  the  line  average  of  solid  density  in  the  fluidized  bed  can  be 
measured  for  each  axial  position  by  recording  the  radiation  attenuation  at 
each  position  and  by  converting  it  to  density  by  equation  (1),  provided  the 
following  assumption  is  mades 


The  fluidizing  particles  and  the  f  iuidizing  gas  together  form  a  substance 
of  density,  j°  s   and  with  a  thickness,  D     (the  inside  diameter  of  the  fluidising 
column),  where  the  absorption  of    #" -rays  by  the  air   (or  fluidizing  gas)  is 
negligible. 

Apparatus  and  Equipment 

A  photograph  of  the  equipment  is  shown  in  Fig*   1  while  a  schematic  diagram 
is  presented  in  Fig.  2.     The  equipment  msg   be  divided  into  two  categories! 
(1),  the  column  and  its  accessories,   and  (2),  the  gamma-ray  source  and  nuclear 
instrumentation  for  detection,  measurement,   and  recording  of  the  transmitted 
gamma-ray  beam. 

The  Column  and  Access or ies.     A  2U- inch-high,   3.97~inch-diameter,   lucite 
column  was  used.     The  column  was  covered  at  the  top  with  a  fine  screen  to 
permit  the  exhaust  of  the  fluidizing  air,  but  to  retain  any  particles  that 
might  have  been  fluidized  as  high  as  the  top  of  the  column.     The  lower  flange 
of  the  column  was  equipped  with  a  pressure  tap  for  measuring  the  pressure  drop 
over  the  bed.     Other  taps  were  located  every  four  inches  up  the  flange  of  the 
column.     The  column  was  bolted  directly  to  the  distributor  assembly. 

Accessories  to  the  column  included  the  distributor  assembly  for  distri- 
buting the  air,  manometers  for  measuring  the  pressure  drops  over  the  column 
and  across  the  distributor,   an  air  supply  to  provide  dry  air  at  a  known 
velocity,  and  a  jack  for  raising  and  lowering  the  column  so  that  different 
heights  in  the  bed  could  be  investigated.     These  accessories  will  be  dis- 
cussed in  this  order. 


Distributor  Assembly.     The  air  was  distributed  by  the  distributor 
assembly  shown  in  Fig.  3.     The  air  first  passed  through  an  8-inch  calming 
section  of  quarter-inch  ceramic  spheres.     Above  the  bed  of  ceramic  spheres, 
the  air  was  distributed  by  a  canvas  filter  cloth  positioned  between  the  lower 
flange  of  the  column  and  a  second  luciie  flange  similar  to  the  lower  flange 
of  the  column. 

Manometers.     A  20-inch  manometer  using  a  manometer  oil  with  specific 
gravity  of  0.818  was  used  for  measuring  the  pressure  drop  over  the  column. 

Air  Supply.     A  compressed-air  line  from  the  Kansas  State  University 
Physical  Plant  provided  air  for  fluidizing.     The  air  was  dried  by  passing  it 
through  a  silica-gel  air  dryer.     The  flow  rate  of  air  was  measured  by  a 
rotameter. 

Jack  Assembly.     The  column-distributor  assembly  was  mounted  on  an  iron 
frame  which  was  rasied  and  lowered  by  means  of  a  hand  jack  making  it  possible 
to  investigate  various  heights  in  the  column. 

Gamma-ray  Source  and  Nuclear  Instrument  at  ion.     Two  gamma-ray  sources 
were  used  to  provide  a  greater  intensity,  thereby  increasing  the  count  rate 

and  improving  the  statistics.     These  two  sources  were  a  five  millicurie 

226  1 "V7 

Ra        source  and  a  twenty-five  millicurie  Cs     Jl   source.     The  two  sources 

were  placed  in  a  5/8-inch  hole  drilled  lengthwise  in  a  2  x  k  x  8-inch  lead 

brick.     The  5/8~inch  hole  served  as  a  collimator  for  the    "#"-radiation  from 

the  two  sources.     For  added  shielding,   other  lead  bricks  were  placed  around 

the  brick  containing  the  source.     The  beam  was  directed  through  the  center  of 

the  column.     On  the  opposite  side  of  the  column  a  7/l6-inch  diameter  collimator 

transmitted  the  beam  to  the  scintillation  probe  which  was  positioned  directly 

behind  the  collimator.     This  collimator  consisted  of  two  bricks.     This  was  done 
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Figure  1.    Photograph  of  column  and  associated  equipment. 
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so  that  the  8-inch  side  of  the  bricks  could  be  used  to  shield  the  probe  from 
stray  radiation.     A  close-up  view  of  the  source,    collimators,   column,   and 
probe   is  shown  in  Fig.   h» 

Instrumentation  for  the  measurement,   detection,   and  the  recording  of  the 
transmitted  gamma-radiation  included  a  scintillation  probe,  power  supply, 
modified  count  rate  meter,  bucking  voltage  supply,   and  recorder „     These 
instruments  will  be  discussed  in  this  order. 

Scintillation  Probe.     The  probe,   a  Model  Dp3,   B-J  Electronics, 
Borg-Warner  Corporation,  was  used  without  modification.     This   instrument 
consisted  of  a  Type  70U,  thallium-activated  sodium  iodide  scintillation 
crystal,   a  Type  6292,   DuMont  photomultiplier  tube,   and  a  one-tube  (6AK5) 
preamplifier.     The  probe  was  operated  at   1200  volts. 

Power  Supply.     A  John  Fluke  Mfg.   Co,   Model  U00BDA  high- voltage  power 
supply  was  used  for  operation  of  the  scintillation  probe. 

Modified  Count  Rate  Meter    (CRM),     The  scintillation  probe  output  was  fed 
into  a  B-J  Electronics  Model  DM1-D  count  rate  meter  that  had  been  slightly 
modified.     The  meter   included  five  time  constants,    0.3,    I,   3,    10,   and  30 
seconds.     Only  the  30-second  time  constant  was  used  for  measurement  of  the 
average  line  density.     The  rate  meter  was   equipped  with  five  scale  multipli- 
cation constants;   10,   30,    100,   300,   and  1000.     300  scale  was  chosen  since  it 
gave  the  highest  sensitivity.     In  choosing  the  300  scale  multiplier  it  became 
necessary  to  modify  the  rate  meter.     With  the  high  activity  used,    it  was 
necessary  to  apply  a  bucking  voltage   in  order  that  the  CRM  would  read 
on-scale,   and  so  that  the  full  range  of  density  fluctuations  could  be 
covered.     The  bucking  voltage  was  added  to  the  input  of  the  vacuum-tube 
voltmeter  circuit  of  the  CRM.     This  bucking  voltage  acted  as  a  negative  bias 
on  the  input  grid  of  the  voltmeter  circuit. 
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Recorder.     This   instrument ,   an  Esterline-Angus  Model  AW  graphic  ammeter 
was  used  without  modification.     The  recorder  was  connected  in  series  with  the 
CRM. 

Preparation  of  Fluidizing  Particles 

The  particles  used  in  this   experiment  were  glass  beads  manufactured  by 
Minnesota  Mining  and  Manufacturing  Co,      In  order  to  minimize  the  error  due  to 
the  size  distribution  effect  when  doing  the   investigation  for  narrow-cut 
particles,   the  purchased  glass  beads  were  reclassified  by  a  sieving  machineo 

The  glass  beads  actually  used  in  the  experiment  were  of  two  sizes i 

(1)  1;0-!|5  U.  S.  mesh:     particle  diameter  ranges  from  0.0138   inches 
to  0.0161).  inches,   the  average  diameter   of  this  size  is  0.015>3 
inches. 

(2)  80-100  U.    S.   mesh:     particle  diameter  ranges  from  0.00U9  to 
0.0070   inches,   the  average  diameter   of  this  size   is  0.006I).2   inches. 

The  particles   of  UO-i.)!?  mesh  were  obtained  by  shaking  about   100  grams   of 
3M  catalog  No.   70  glass  beads   on  the  sieving  machine  for  three  minutes,   those 
particles  passing  through  kO  mesh  screen  and  retained  on  the  h$  mesh  screen 
were  collected. 

The  80-100  mesh  particles  were  obtained  from  the  catalog  No.    110  glass 
beads  by  using  the  same  technique  as  for  bP-U5  mesh  particles.     Those  passing 
through  80  mesh  and  retained  on  100  mesh  particles  were  collected. 

Some  additional  data  were  obtained  with  60-70  mesh  particles   in  order  to 
confirm  some  of  the  observations  and  correlations  resulting  from  the  use  of 
the  other  two  particle  sizes. 
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Preliminary  Calibrations 

Cotmt  Rate  Meter  vs.  Solid  Dens ity  in  the  Fluidized  Bed.     For  the  prelimi- 
nary calibration,  spacers  made  from  very  thin  aluminum  sheets  were  used,  which 
had  inside  air-gap  widths  of  0.25,  0.50,  0.75,   1.00,   1.25,   1.5>0  &  2.00  inches, 
respectively. 

The  calibration  was  done  by  inserting  perpendicularly  one  spacer  at  a 
time  into  the  densely  packed  bed  of  glass  beads  of  density  equal  to  1.502 
grams/cu-cm,   and  then  recording  the  radiation  attenuation.     This  process  was 
also  repeated  by  packing  the  glass  beads  inside  the  spacers,  and  leaving  the 
area  surrounding  the  spacers  empty.     Then  the  density  corresponding  to  a 
particular  spacer  could  be  calculated  with  the  following  two  relations? 

(a)  (Density  of  packed  bed)  x  (D    -  D  )/D  ,  when  the  glass  beads 

c    s   c 

were  outside  the  spacers. 

(b)  (Density  of  packed  bed)  x  (Ds/Dc),  when  the  glass  beads  were 
inside  the  spacers. 

These  two  linear  relationships  can  be  proved  by  the  combined  Lamberts  equation: 

m(i/i0)  -  <Mjcfc  d0  ♦  Ctf«)jufiu  \i  *  U/.WV  £        (2) 

Where  the  subscript  »c"  denotes  the  radiation  attenuation  due  to    y-rayS 
absorption  by  the  walls  of  column^  subscript  RA1",  by  the  aluminum  sheets  and 
subscript  "p",  by  the  packed  particles  inside  spacers  or   in  the  annular  space 
between  column  walls  and  the  spacer  sheets.     Care  was  taken  that  the  attenua- 
tion due  to  column  wall,  aluminum  sheets  and  all  other  possible  absorbing 
media  stayed  constant  through  the  calibration  (Fig.  k  and  5)»     Therefore, 
Equation  (2)  can  be  simplified  ass 
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in  (I/IQ)  p  constant  *  (yUm)pfpa  ^  (3) 

where  (y£/m)p  P-na  *s  a*so  constant  throughout  the  calibration,  since  the  same 
glass  beads  were  used,  and  the  density  of  solid  glass  remained  constant.  Then 
the  change  in  ln(l/l0)   is  a  linear  function  of  D_ 

In  dealing  with  the  density  of  a  fluidizing  state,  the  radiation  attenua- 
tion reading  corresponding  to  one  of  the  spacer  calibrated  values  should  have 
a  density  corresponding  to  the  mass  of   f>Da  D  P^     particles  distributed 
homogeneously  in  the  volume,  DcAr,  where, 

A_,  ■  projected  area  of  y"-beam  path. 
Therefore,  the  density  of  the  fluidized  bed  at  a  certain  attenuation  reading  iss 

1°    -Jpa  Dp  ^r     J  pa  \  „  N 

where,  Dp     s  Dr  -  Ds  ,  when  the  beads  were  outside  the  spacer . 
Dp*   ■  Ds  ,  when  the  beads  were  inside  the  spacer. 

As  described  in  the  section  of  Apparatus  and  Equipment,   the  radiation 
attenuations  were  measured  and  recorded  on  the  graphic  ammeter.     In  using  the 
300  scale  multiplier  of  the  count  rate  meter,    it  was  necessary  to  calibrate  the 
density  vs»  recorder  reading  in  the  manner  mentioned  above  for  several  different 
bucking  voltages.     This  was  done  so  that  the  full  range  of  densities,   i.e.,  from 
empty  column  to  packed  bed,   could  be  covered. 

A  complete  correlation  chart  of  recorder  reading  vs.  density  of  bed  is 
given  in  Table  I.     In  the  tables 

300  I  scale  is  for  B.  V.   of  8170—8300  v 

300  la  scale  is  for  B*  V,  of  9870—10010  v 

300  II  scale  is  for  B.  V.   of  11300— 11U50  v 


/♦ 


Figure  4.    Bed  arrangement  with  spacer  out. 
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Figure  5.    Calibration  technique  with  spacer  inserted. 
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300  Ha  scale  is  for  B.  V.   of  12100—12220  v 

300  III  scale  is  for  B.  V.   of  13100—13330  v 

Fig.  6  and  Fig.  7  are  two  typical  tracings  of  the  CRM  recorder  on  300  I 
scale  and  300  III  scale  respectively. 

Rotameter  Calibration.     The  linear  flow  velocity  of  the  air  was  measured 
by  a  rotameter  of  Brooks  Rotameter  Co,  R-9M-25-2.     Two  floats  and  one  pressure 
gauge  were  used  to  cover  the  range  of  flow  rates  used  in  this  experiment.     The 
aluminum  float  was  used  to  measure  the  low  flow  rates,  and  the  steel  float  was 
used  to  measure  high  flow  rates.     The  rotameter  was  further  calibrated  by  an 
Anemotherm  (Wo.   1005  Mechanical  Engineering.,  K.S.U.)  for  the  high  flow  rates. 
For  calibration  in  the  low  range,,   a  wet  test  meter   (WMT  JLj...  Dept.   of  Chemical 
Engineering)  was  used.     Five  pressures  of  $,   10,   l£,   20^   2$  psig  were  used  as 
parameters.     The  calibration  curve  is  presented  in  Fig.   8. 

Procedure 

The  instruments  were  turned  on  at  least  it  hours,  usually  overnight, 
before  a  run  was  ma.de.     The  calibration  was  checked  to  make  sure  that  the 
instrumentations  were  stabilized.     The  column  was  charged  with  the  desired 
particles.     Pressure  drop  data  were  taken,,    increasing  the  air  velocity  from 
below  the  minimum  fluidizing  velocity  to  about  90  ft/min.     The  air  velocity 
was  then  decreased  slowly  until  the  air  velocity  was  zero.     The  static  bed 
height  was  then  read.     The  air  was  then  set  at  the  desired  velocity.     The  data 
were  taken  with  300  multiplier  on  30-second  time  constants  and  with  the  bucking 
voltage  appropriate  for  the  particular  density  range.     Data  were  generally 
taken  at  every  half  inch  in  the  bed,   starting  1  inch  above  the  distributor. 


2U 


Near  the  top  of  the  fluidized  bed,  data  were  taken  every  quarter  inch.  The 
attenuation  readings  from  the  recorder  chart  were  converted  to  average  line 
density  by  the  calibration  table  (Table  I),  plotted  on  rectangular  graphs  to 
obtain  the  density  profile,  and  then  the  consistency  test  performed. 

Consistency  Test 

To  prove  that  the  radioactive  technique  used  for  measuring  the  axial 
solid  density  was  satisfactory  and  reproducible,  the  consistency  test  for 
several  sets  of  experimental  data  was  performed,,  The  test  was  based  on  the 
principle  of  the  conservation  of  mass,  i.e.,  the  actual  mass  of  particles 
placed  in  the  column  should  be  equal  to  the  integrated  area  of  the  solid- 
density  profile. 

For  most  runs  the  results  obtained  hy   this  consistency  test  were  within 
the  allowable  error.  The  procedure  and  the  calculation  method  of  this  consis- 
tency test  are  as  follows: 

The  area  under  the  curve  shown  in  Fig.  9  was  graphically  integrated  by 

means  of  a  planimeter.  The  total  mass  of  particles  determined  in  this  manner 

is: 

M*  «  (integrated  area)  x  (cross  sectional  area  of  column)  x 
(unit  conversion  constant) 

If  the  actual  mass  of  particles  in  the  column  is  M,  then  the  deviational  error 

%   error  - 5  x  10C$ 

-M 

Table  2  summarizes  the  results  of  this  consistency  test.  The  positive  or 
negative  error  for  each  run  is  listed. 
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The  factors  which  caused  the  error  of  this  consistency  test  may  be  attri- 
buted to  j 

(a)  The  assumption  that  the  average  line  density  of  solid  is  uniform 
across  the  cross  sectional  area  at  each  axial  position  is  not  exactly  valid 
since  the  density  near  the  wall  of  column  could  be  different  from  that  in  the 
middle  of  the  column. 

(b)  There  was  an  error  due  to  radioactive  decay  or  noise  and  other  effects 
on  the  detecting  power  of  the  scintillation  probe  for  "^T-ray  attenuation. 

(c)  When  fluidizing,  the  particles  in  the  bed  were  actually  in  a  state  of 
dynamic  equilibrium.  The  line  density  detected  is  an  average  within  the  time 
constant  used  for  the  scintillation  counting  probe.  Because  of  the  complexity 
of  the  inter-particle  disturbance,  the  speed  of  the  particles  at  different 
locations  in  a  certain  layer  could  not  be  the  same,,  but  the  time  constant  usea 
for  the  scintillation  counting  probe  was  the  same  throughout.  Therefore  it 
could  hardly  be  possible  to  take  an  exact  average  solid  density  which  was 
consistent  at  any  cross-section  of  the  bed» 

(d)  The  graphical  integration  must  have  introduced  some  error. 

Fig*  10  shows  the  positive  or  negative  percent  error  of  this  consistency 
test  as  a  function  of  air  velocity.  The  negative  error  consistently  increased 
as  the  flow  rate  of  air  increased.  At  low  flow  rate  of  air,,  the  error  was 
almost  negligible. 

The  fact  that  a  larger  error  occurred  at  higher  air  velocity  is  logical 
since  the  particles  were  moving  more  vigorously  and  randomly  in  the  bed.  It 
might  also  be  due  to  elutriation  of  small  particles  from  the  open  top  of  the 
fluidizing  column. 
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CHARACTERISTIC  AXIAL  DENSITY  DISTRIBUTION 
OF  GAS-SOLID,  AGGREGATIVE  FLUIDIZED  BEDS 


General  Observations 


Typical  characteristic  axial  solid  density  profiles  of  a  gas-solid 
fluidized  bed  are  shown  in  Fig.  11  and  Fig*  12.  Similar  profiles  were 
obtained  for  all  experimental  runs  for  different  static  bed  heights  and  under 
different  air  velocities.  Even  the  beds  mixed  particles  yielded  similar  pro- 
files. Two  distinct  zones  can  be  found  in  this  typical  distribution  curve: 
in  the  range  AB  of  Fig.  11,  the  solid  density  stays  fairly  constant.  This 
is  called  the  constant  density  zone.  After  point  B  the  solid  density  decreases 
sharply  as  the  distance  from  the  bottom  of  the  bed  increases.  This  portion  of 
the  profile  is  called  the  falling  density  zone. 

These  typical  density  profiles  show  qualitatively  the  properties  of  gas- 
solid  fluidized  beds  under  the  effect  of  several  operational  variables:  — 

Dotson  (5)  studied  the  effects  of  five  primary  variables  on  the  uni- 
formity of  fluidization  in  a  U-inch  diameter  column  by   a  statistically 
designed  experiment.  He  measured  the  density  fluctuations  inside  the  bed  by 
a  capacitance  method  and  concluded  that  among  the  five  variables,  i.e., 
distributor,  bed  height,  particle  size,  probe  height  and  gas  velocity,  gas 
velocity  had  by  far  the  greatest  effect  on  the  uniformity,  with  uniformity 
generally  decreasing  with  an  increasing  gas  velocity. 

In  the  present  study  with  the  ^T-ray  technique,  the  capacitance  probe, 
which  did  effect  the  flow  pattern  in  a  fluidized  bed,  was  not  present.  Thus 
its  effect  was  eliminated.  Only  one  type  of  gas  distributor  (canvas  filter 
cloth)  was  used  in  all  the  experimental  runs.  Therefore,  only  three 
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operational  variables  were  studied.  The  effect  of  these  three  operational 
variables  on  the  axial  solid  density  distribution  -will  be  discussed. 

Effect  of  Air  Velocity 

Figo  11  shows  the  effect  of  air  velocity  upon  the  fluidized  bulk  density 
profiles.  It  is  shown  that  the  density  in  the  constant  density  zone  decreased 
with  increasing  air  velocity,  but  the  slope  of  the  falling  density  zone  becomes 
smaller  as  the  air  velocity  increases.  This  Xiras  expected  since  more  air  was 
passing  through  the  bed  at  higher  air  velocity  and,  therefore,  the  bed  was 
expanded. 

Fig*  11  also  indicates  that  this  two-zone  profile  becomes  less  apparent 
as  the  air  flow  rate  increases.  At  flow  rate  over  90  ft/min.,  the  line  average 
density  falls  from  the  bottom  of  the  bed  to  the  top  layer  of  the  fluidizing  bed. 
In  other  words,  when  the  air  velocity  is  low  the  constant  density  zone  clearly 
exists  and  the  bed  is  closer  to  the  ideally  expanded  fixed  bed.  If  the  air 
velocity  is  high,  the  air  which  carries  the  particles  becomes  more  turbulent 
and  the  particles  are  in  more  vigorous  and  randomized  movement.  Under  these 
conditions  there  is  no  section  in  the  whole  bed  that  has  a  uniform  solid 
density. 

Effect  of  Static  Bed  Height 

Figs.  13  (a-d)  are  the  dimensionless  axial  density  profiles,  in  which 
f/f         was  plotted  against  h/h   .  It  appears  from  these  graphs  that  the 
static  bed  effect  on  the  density  profile  is  relatively  small.  The  dimension- 
less  density  profiles  for  different  static  bed  heights  fluidizing  at  the  same 
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air  velocity  can  almost  he  represented  hy  one  curve.     This   is  especially  true 
at  low  fluidizing  air  velocities.     The  deviations  between  such  curves  becomes 
greater  as  the  air  velocity  increases.     The  deviation  from  the  common  dimen- 
sionless  density  profile  becomes  more  significant  for  the  individual  profile 
corresponding  to  the  smallest  bed  height.     This  is  shown  in  Figs.   13   (c)  and 
13  (d).     This  result  may  be  explained  as  follo-nis: 

The  ratio  of  the  static  bed  height  to  the  column  diameter  instead  of 
the  static  bed  height  itself  is  probably  the  significant  variable  affecting 
the  bed  expansion  and  the  density  profile.     This  is  the  proper  combination 
of  entrance  and  wall  effects.     When  the  bed  is  shallow  and  the  ratio  is  small, 
the  air  passing  through  could  hardly  have  time  to  form  well-shaped  bubbles 
before  it  leaves  the  fluidizing  bed.     The  two  phases  of  fluid izat ion  for  the 
shallow  bed  are  therefore  not  as  pronounced  as  in  the  case  of  deep  static  bed. 

It  may  be  concluded  that  at  higher  air  velocity,   the  effect  of  the  static 
bed  height  on  the  densitjr  profile  becomes  more  noticeable  when  the  ratio  of 
static  bed  height  to  the  column  diameter  is  less  than  one  (the  diameter1  of 
fluidizing  column  used  was  3.9765  inch).     If  the  ratio  is  greater  than  one, 
the  effect  of  the  static  bed  height  on  the  dimensionless  density  profile  is 
negligible.     This  is  consistent  with  the  observation  made  by  Gomezplata 
et.al.,(8).     They  studied  the  entrance  effect  on  the  conversion  of  the 
catalytic  cracking  of  cumene.     Their  results  are  shown  in  Fig.    Ik  to  compare 
with  and  to  substantiate  the  present  observation. 
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Effect  of  Particle  Size  and  Bed  Composition 

The  effect  of  particle  size  and  bed  composition  could  not  be  determined 
either  qualitatively  or  quantitatively  without  defining  some  index  for 
aggregative  fluidization.  Stewart  (25)  has  discussed  the  particle  size 
effect  with  a.  defined  ISU  index.  His  results  are  consistent  with  the 
findings  of  Dotson  (5)  and  Baumgarten  (3)» 

More  details  and  quantitative  discussion  of  the  effects  of  operational 
variables  on  aggregative  fluidization  will  be  found  in  later  sections. 

ANALYSIS  CF  DATA  AND  CORRELATIONS 

As  mentioned  in  the  previous  section,  the  axial  density  profiles  may  be 
divided  into  a  constant  density  zone  and  a  falling  density  zone.  The  corre- 
lation of  the  density  profile  as  a  function  of  the  operational  variables  will 
be  obtained  for  the  two  zones  separately.  It  was  usually  considered  (1)  that 
the  constant  density  zone  extended  from  the  bottom  of  the  bed  to  the  position 
corresponding  to  the  static  bed  height,  and  the  falling  density  zone  existed 
above  the  static  bed  height. 


Bed  Density  Correlation  in  the 
Constant-Density  Zone 


The  density  in  the  constant-density  zone  is  fairly  constant,  thought  a 
considerable  density  fluctuation  may  be  observed  from  the  density  profile  at 
high  air  velocities.  From  the  qualitative  observation  made  in  the  previous 
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section,  we  know  that  air  velocity  is  the  major  factor  affecting  the  shape  of 
the  density  profile,  whereas  the  static  bed  height  has  little  effect.  Thus 
the  density  in  the  constant- density  zone  can  be  correlated  as  a  function  of 
air  velocity  and  particle  diameter. 

Let  Pft    represent  the  average  density  in  the  constant  density  zone 
(where  the  region  ranges  from  the  bottom  of  the  column  to  the  position 
corresponding  to  the  static  bed  height,  1^)3  Ms  ,  the  weight  of  particles 
fluidizing  below  L^;  and  A,  the  cross- sectional  area  of  the  fluidizing 
column,  thent 

?< "  It  (5) 

The  particles  blown  off  and  suspended  above  L^,  are  then 

M  -  HSl  +  M^  (6) 

where  M  is  the  total  weight  of  particles  charged  in  the  column.  It  was 
mentioned  by  Bakker  et.al.  (2)  that  the  volume  of  the  particles  suspended  in 
the  space  above  JLjnf,  i.e.  Vols  ,  corresponds  to  the  product  of  the  volume  of 
the  gas  bubbles  appearing  below  L-,  and  the  concentration  in  the  prefluidized 
state.  Thus, 

VolSi  -  (l-fnf)  Volb  (7) 

If  we  consider  that  A'  is  the  part  of  area  A  which  is  occupied  by  the  dense 
phase  emulsion,  then  (A  -  A1)  is  the  area  crossed  by  the  bubbles.  The  super- 
ficial velocity  of  bubbles,  as  indicated  hy  Baumgarten  et.al.  (3),  is  a 
function  of  bubble  size  and  is  different  at  each  layer  of  the  fluidized  bed. 
He  observed  that  the  size  of  a  bubble  grows  as  it  moves  away  from  the 


*  £mf>  *s  ^e  vo^!^  space  of  bed  at  minimum  fluidizing  state. 


hh 


distributor.  But  to  avoid  this  complexity,  we  can  assume  that  there  ic  a  mean 
bubble  velocity,  VY,,  which  is  constant  throughout  the  section  below  L^.  With 
t h i s  as sumpt ion: 

G  -  %  *   G^  (8) 

where  G  is  the  volumetric  air  velocity,  Gb  is  the  volumetric  bubble  velocity, 
and  Gj^,  the  volumetric  velocity  of  minimum  fluidization. 
Therefore,  from  equation  (8)s 

VA  =  Vb  (A  -  A»)  +  Vmf  A»  (8») 

Rearranging  equation  (8')  we  get; 

A  (Vb  -  V) 
Vb  "  Vaf 

Let  tb  be  the  residence  time  of  bubbles  in  the  bed  below  L-.  Then: 

W 
tb  -  —  (10) 

vb 
The  volume  occupied  by  bubbles  below  L^  ist 

Lmf 
Volb  ^  Gb  tb  -  Vb  (A- AOj-  (11) 

Substituting  equation  (9)  into  Equation  (11)  gives: 

A(?b  -  V) 


Volb  "  Lmf 


A  - 


I     (Vb  -  Vaf ) 


_  ^mf  (V  -  Vmf ) 

%  ~  V  (12) 

From  Equations  (7)  and  (12),  the  volume  of  the  solid  particles  above  L-  is 
obtained. 


h$ 


Al^f  (V  -  V_,f ) 

"V'T  v     (1-^}  (13) 

Vb  "  vmf 

then, 

AL  ,.  (V  -  V  -) 

'        x   J  vb  ~  v^ 

and 

2  vb  -  Vmf 

Substitution  of  Equation  (1J?)  to  Equation  (5)   leads  to: 

fA  -  *L  -  JL_  .    /sd-g-f)  (v  -  V)  (16) 

^mf       ALmf  (Vb  -  Vmf) 

since, 

1        £  mf  ^"p— 

then, 

"    ?  mf  -  C<V  "  Vmf)  (17) 


where     C  ■»  -i-i°£. 


<Vb  "  \f ) 

Since  the  average  bubble  velocity  through  the  bed  below  L^.  has  been  assumed 
constant,   C  in  equation  (17)   is  a  constant  and  the  average  density  in  the 
constant  density  zone  has  a  linear  relationship  with  the  superficial  air 
velocity. 


he 


The  data  plotted  on  Fig.   lf>  (a)  and  Fig.   15  (b)  confirm  this  linear 
relationship  between    P  .  and  (V  -  V  - )  fairly  well  with  the  static  bed 
height  showing  a  negligible  effect. 

Using  the  least  squares  method,  the  intercepts  and  slopes  of  the  straight 
lines  from  the  experimental  data  were  obtained.     The  empirical  equations 
which  relate    P  d  with  (V  -  Vmf)  for  hP-k5  raesh  and  80-100  mesh     3M  glass 
beads  are  as  follows: 

(i)       For  J4.O-U5  mesh  glass  beads: 

f d  =  1.371  -  0.00536  (v  -  v .  h0  mesh) 

(ii)     For  80-100  mesh  glass  beads j 

fd  ■  i'26?  -  °-00^  <V  -  V   80  nesh' 
where     P  3  is  in  gm/cu-cm 

V,  Vj-f  are  in  ft/min 

^mf     hO  mesh  *s  ^ie  m^-n^mm  velocity  of  fluidization  of  iiO-U5 
mesh  glass  beads  with  air,  which  is  2U. 3  ft/min 

V_f     qq  jnggh  is  the  minimum  velocity  of  fluidization  of  80-100 
'  mesh  glass  beads  with  air,  which  is  I4..80  ft/min 


Bed  Density  Correlation  in  the 
Falling  Density  Zone 


The  fluidized  bulk  density  in  this  zone  decreases  as  the  height  from 
distributor  increases.  In  all  cases,  as  has  been  shown  in  Fig.  13,  the 
density  distribution  in  this  falling  density  zone  may  be  represented  by  a 
probability  function.  For  the  experimental  runs  at  low  air  velocity,  the 
density  distribution  in  this  falling  zone  can  be  fitted  almost  exactly  by 
the  one-tailed  normal  distribution  function,    if  we  consider  that  the 
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falling  density  zone  starts  at  the  location  corresponding  to  the  static  bed 
height.  But  for  the  falling  density  distribution  obtained  from  the  bed 
fluidizing  at  higher  air  velocity,  the  profile  becomes  rather  skewed.  We 
can  hardly  find  a  simple  form  of  probability  equation  to  fit  this  profile. 
Consequently,  it  is  not  an  easy  task  to  find  a  unique  and  general  form  of 
probability  equation  capable  of  describing  the  falling  density  profiles  for 
all  ranges  of  air  velocity. 

Thus,  it  becomes  necessary  to  search  for  an  empirical  equation  which, 
with  some  coefficients  determined  from  the  experimental  data,  will  describe 
quantitatively  the  state  of  aggregative  fluidization.  For  instance,  the 
varying  slope  of  the  falling  density  profile  could  be  a  simple  index  for  this 
purpose. 

It  was  found  that  a  function  in  the  form: 

f  //pa  "  *»  "*  C  B°  (hA»*  -  1)2  3  (18) 

would  correlate  fairly  well  the  falling  density  profiles.  At  a  height  corres- 
ponding to  L  .  (the  position  at  prefluidized  state)  the  fluidized  bulk  density 

mf 

is  P  3,  which  can  be  obtained  from  the  axial  density  profile.  In  most  cases 
it  was  the  same  as  9  d,  the  average  density  in  the  constant  density  zone. 
Values  of  Pd  have  been  plotted  as  function  of  air  velocity  and  particle  size 
in  Fig.  l£.  The  value  of  A0  in  equation  (18)  is  then  equal  to  J°d/ppa> 
which  can  be  calculated  for  each  condition. 

The  problem  now  is  to  solve  for  B0,  which  indicates  quantitatively  the 
slope  of  the  correlated  density  profile  for  the  falling  density  zone.  It  may 
be  considered  as  an  index  to  describe  the  fluidizing  states  under  various 
operating  conditions. 
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By  taking  logarithms  of  both  sides  of  equation  (18),  we  get  the  following 
linear  equation; 

2 


log  (f/P)   =  log  (JVf  pa>  ♦  B  (h/Lmf  -  1)' 


pa 


=  A  +  B  (h/L  .  -  l)2  (19) 


It  should  be  noted  that  equation  (19)  is  only  an  empirical  equation  to  fit 
the  experimental  falling  density  profiles.  It  is  not  fitted  equally  well  by 
the  experimental  data  obtained  under  different  fluidizing  conditions,,  The 
general  trend  is  shown  in  Fig.  16. 

For  all  sizes  of  glass  beads  fluidizing  at  lower  air  velocities,  the 
falling  density  profiles  were  fitted  better  by  the  linear  relation  of  equation 
(19);  xtfhile  for  those  fluidizing  at  higher  air  velocities,  the  linear  relation- 
ship was  not  the  best  fit  to  the  profile.  Fig.  16  is  the  falling  density 
profiles  transformed  into  the  form  of  equation  (19) «  Even  though  the  linear 
correlation  of  the  dimensionless  density  to  the  dimensionless  height  was  not 
as  good  for  the  runs  fluidizing  at  higher  air  velocities  as  that  at  lower 
velocities,  it  is  possible  to  obtain  a  best  straight  line  by  using  the  method 
of  least  squares.  Thus  the  slope,  B,  and  the  intercept,  A  =  log  (  PVP  )  can 
be  found  for  all  sets  of  experimental  data  fluidizing  at  various  conditions. 
They  are  listed  in  Table  3  (Appendix). 

The  index  B  was  plotted  in  Fig.  17  (a),  17  (b)  and  17  (c)  with  the  air 
velocity  variable  for  three  different  sizes  of  glass  beads.  This  index 
converges  to  a  single  value  as  the  air  velocity  increases.  It  indicates  that 
at  high  air  velocity,  the  particles  carried  up  by  the  air  bubbles  are  uniformly 
dispersed  in  the  falling  density  zone.  Wo  matter  what  the  original  static  bed 
height,  the  value  of  index  B  is  the  same  if  fluidizing  at  high  air  velocity. 
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Fig.    18  further  verifies  that  the  effect  of  the  static  bed  height   on  aggre- 
gative fluidization  is  pronounced  only  at   low  air  velocity  close  to  the 
minimum  velocity  of  fluidization.      In  Fig,    19,   the  plot   of   index  B  vs.    the 
particle  size  shows  that  the  particle  size  affects  the  state   of  aggregative 
fluidization  only  at   lower  air  velocities.      In  other  words,   the  particles  are 
distributed  in  a  similar  manner   ever  the  falling  density  zone   (at   high  air 
velocity),   no  matter  what  size  of  particles  the  bed  is  composed  of. 

The  value  of  A,  which  is  actually  the  ratio  of  fluidized  bulk  density 
at  the  position  of  the  static  bed  height  to  the  bulk  density  of  the  same 
particles   in  their  packed  state,  was  correlated  in  a  similar  way  as  that 
done  for   index  B.     Fig.   20  shows  that  the  index  A  is  varying  with  air  velo- 
city for  all  static  bed  heights.     In  all  cases   (Fig.   20,   a-c)  the  effect  of 
static  bed  height  on  the   index  A  is  significant  when  Lmf/D   <C  1.     When  this 
ratio   is  greater  than  one,   the  fluidized  bed  expansion  was   essentially  the 
same  for  all  static  bed  heights  fluidizing  at  the  same  air  velocity.     For 
the  80-100  mesh  glass  beads,   the   index  A  correlation  with  air  velocity 
converges  to  one  curve  for  all  static  bed  height,    even  when  L  f/D     is   less 
than  one.     This   leads  to  the  conclusion  that  when  the  particles  size  is  small, 
the  entrance  effect  does  not  affect  the   index  A  (Fig.    20,    c). 

The  separate  effect,  of  particle  size  on  the   index  A  is  mostly  due  to  the 
relative  magnitude  of  the  bulk  density  of  different  particles   in  the  packed 
state  (    P-)a)»     The  size  effect   is  practically  the  same  at  all  air-velocity 
levels  as  has  been  shown  in  Fig.    21. 

The  effect   of  bed  composition  on  the  values   of   index  A  and  index  B   is 
shown  in  Fig.   22.     The  beds  were  originally  at  a.  static  height   of  6.5  inch  and 
were  composed  of  various  fractions   of  hO-hS  mesh  and  80-100  mesh  particles. 
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They  were  fluidized  at  an  air  velocity  of  60  ft/min.     The  axial  density  pro- 
files were  obtained  and  index  A  and  index  B  calculated  for  every  bed  composi- 
tiono     Bed  compositions  of  0,   2,   £,    10,   £0,    90,   95,   98  and  100  percent  by 
weight  of  L|.0~1|£  mesh  glass  beads  were  used.     The  results  show  that  the  mixed 
particle  size  affects   index  A  and   index  B  significantly  only  when  the  composi- 
tion  is  below  10%  or  above  90%  by  weight  of  h0-h5  mesh  glass  beads . 

As  the  state  of  aggregative  fluidized  beds   is  affected  by  the  three  major 
variables,    i.e.,   superficial  fluid  velocity,  particle  size  and  the  ratio  of 
static  bed  height  to  the  diameter  of  column:;    index  B  as  well  as   index  A 
(which  is  used  to  describe  the  state  of  aggregative  fluidization)   should  be 
correlated  in  a  single  empirical  equation.     By  applying  the  Pi-theorem  in 
the  dimensional  analysis,  we  obtains 

I  A  I  -  Ol  '  (  ^P~)  a  <  VUb  <WA>c)C  <2°) 

where,   the  values   of    0(  ',   a,   b,   cj    0    ',   a',   b'   and  c'   can  be  obtained:  from 

the  experimental  data. 

The  equations  for  the  series   of  data  obtained  ares 

-S/DnV  ft\0.978  -1.265  O.3U69 

I  A  I  -     1.396  x  10  \  JLJJL)  (Dp/Dc)  (Lmf/Dc)  ^ 

9/DnV  f/=\-2.09li  3.0U5  0.3ii7U 

J  B   J   -     5.U92  x  10  (^TJ  W  (Lmf/Dc)  ^ 

Equations  (22)  and  (23)  are  plotted  in  Fig.  23  and  Fig.  2h  in  order  to 
compare  the  calculated  values  for  |  B  |  and  |  A  I  with  those  experimentally 
determined  (Table  h  and  Table  5). 
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Table  1*.     Comparison  of  calculated     |B  |    with  the  experimental   data. 


1*0-1*5  mesh  g.l 

.ass  beads              i 

!               80-100  mesh  glass  beads 

Experimental       : 

Calculated       : 

Experimental        :       Calculated 

16.789 

12.521* 

25.695                         18.910 

3.1*38 

2.932 

2.925                           i*.  382 

1.391 

1 .  26U 

1.572                           1.391 

0.813 

0.830 

1.710                           2.231 

21.U71 

16.009 

5.602                           6.298 

3.151 

3.778 

2,^50                     2.1*96 

1.61*5 

1.635 

2.035                           1.560 

1.133 

1,071; 

6.192                          6.660 

36.981; 

18.51*8 

60-70  mesh  glass  beads 

5.20U 

ii.35'8 

2.6^0                          2.637 

1.213 

1.861* 

2.11*2                          1.651 

0.815 

1.253 

2.11*1                          5.929 

10.395 

13.280 

1.1*06                          1.388 

3.328 

3.102 

1.135                          0.59li 

1.033 

1.261* 

1*.766                          7.61*2 

2U.010 

17.1*1*1* 

2.185                          1.762 

3.151 

U.051 

1.238                          0.762 

1.150 

1.710 

Table  5.-   Comparison  of  calculated    |  A  |    with  the  experimental  data. 


1*0-1*5  mesh  gla 

ss  beads 

1                80-100  mesh 
l        Experimental 

.  glass  beads 

Experimental        s 

Calculated 

Calculated 

0.0535 

0.0523 

0.0793 

0,0701* 

0.1052 

0.1030 

0.1667 

0.1386 

0.1722 

0.1526 

0.207U 

0.2061 

0.2072 

0.1857 

0.1055 

0.0868 

0.0691 

0.0705 

0.1663 

0.1723 

0.1.1*52 

0.1381* 

0.1993 

0.251*8 

0.21*23 

0.201*7 

0.0983 

0.1010 

0.2789 

0.21*91 

0.1621* 

0.1990 

0.0677 

0.0800 

0.1525 

0.1571* 

60-70  mesh 

glass  beads 

0.2617 

0.23U1 

0.111*2 

0.1625 

0.2912 

0.2818 

0.0991 

0.1301* 

0,01*63 

0.0551* 

0.1318 

0.0899 

0.0825 

0.1092 

0.li*20 

0.1386 

0.1575 

0.1661 

0.1716 

0.1726 

0.0689 

0.0702 

0.1029 

0.0961 

0.1223 

0.1388 

0.1629 

0.11*82 

0.2151 

0.2072 

0.1987 

0.181*1* 

0.0802 

0.0806 

0.1? kh 

0.1596 
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We  can  now  calculate  a  complete  axial  density  profile  for  the  air-particle 
fluidizing  system  with  Fig.    15  and  equations  for    |  A  |     and    |  B  |     .     It   is 
illustrated  in  the  following   illustration: 
Illustration: 

In  a  fluidization  test,    glass  beads  of  IjO-lifJ  mesh  were  used  (Run.     PJ4.O-O6), 
experiment  was  run  at  room  temperature  and  with  air  at   1  atmospheric  pressure 
in  a  circular  column  of  h  inch  diameter.     The  bed  had  a  height  of  6.20  inch  at 
incipient  fluidization.     Estimate  the  fluidized  bulk  density  profile,   for  the 
bed  fluidizing  at  1  ft/sec.   of  superficial  air  velocity.     Compare  it  with  the 
experimental  data. 
Solutions 

Properties  of  air  at  the  fluidizing  condition? 
J°F  *  0.0757  lb/cu-ft 
yd      ■  1.23  x  10"^     lb/ft-sec.  (11) 

The  air  velocity  at   incipient  fluidization  can  be  calculated: 

Vmf  "  22,3^  *Vmin. 
Properties  of  particles: 
D     ■  0.0153  inch 

J3  pa  "*  U*02     Wcu-cm 
Constant  density  zone; 

V  a  1  ft/sec.  *  60  ft/min. 

V  -  V^  *  37.65  ft/min. 

from  Fig.   15  (a),       Pd  ■  1.168  gm/cu-cm. 
Falling  density  zones 

V/f    Dp       (1)(0.0757)(0.0153)       _  M 

NRo  a 77 «  p =  7»o5 

Ke  S*  (1.23  x  l(f5)(12) 
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0.978  n    -1.265  t,  0.31*69  0-CT8o.01ft  -1.265  ,200.31|69 

(KdJ         W  (5T}  (?,85)        (~T")  hf0 


Re' 


*  (7.U9)(1.120)(1.16U)  «  9.763  x  10 
From  Fig.    23, 
|  A   |     -  0.139 

A     -  log  (Pd/J>pa)  ■  -0.139 
-2.09U  L3.0WL,  0.31*74  R.  ~2.091*  0>0153  3.01*5  6.20n0.3U71* 


k  ' 


o  -8 

■  (1.31*  x  10     )(1*. 335  x  10     ) ( 1 .  167 ) 

-10 
*  6.765  x  1.0 

From  Fig.   21*, 

B  |   =  3.65  B  *  -3.65 

Thus,   the  density  in  the  falling  density  zone  can  be  calculated  by  the  follow- 
ing equation; 


l0g  {f/fw}  '  A  +  B  (h/Lmf  ■  1} 


i  .e, 


log  (  J^/1.502)  =  -0.139  -  3.65  (h/6.2  -  1) 
The  calculated  density  and  the  experimental  values  for  the  falling  density  zone 
are  listed  as  follows: 


Height  from     : 
distributor 
h,    in           : 

Experimental 

Calculated 

f//pa     ! 

gm/cu~cm 

1    f/f-pa     '          gm/cu-cm 

7.00 

8.00 

9.00 

10.00 

11.00 

0.611* 
0.31*1* 

0.153 
0.039 
0.020 

0.922 
0.516 
0.229 
0.058 
0.030 

0.633                  0.951 
0.358                  0.538 
0.132                  0.200 
0.031                   0.01*7 
0.005                  0.007 

The  experimental  data  are  plotted  in  Fig.  ?5  with  the  density  profile  cal 
culated  by  the  obtained  equation. 
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Determination  of  the  Fluidized  Bed  Height 

The  behavior  of  gas-solid  fluidized  beds  observed  in  this  work  deviated 
considerably  from  the  behavior  of  an  ideal  expansion  of  a  packed  bed.  The 
particles  were  not  as  'uniformly  suspended  in  the  bed  as  is  usually  observed 
in  most  of  the  liquid  fluidized  beds.  Even  the  two  phase  model  for   gas-solid 
fluidization  is  only  a  very  approximate  one.  Actually,  there  is  no  distinct 
phase  which  can  be  defined  as  a  bubble  phase  consisting  of  purely  gases,  nor 
is  there  a  dense  phase  consisting  of  uniform  emulsion  of  particles.  The 
bubbles,  formed  by  the  gas  and  held  in  the  pore  spaces  between  solids,  grow 
in  size  as  they  rise  from  gas  distributor  toward  the  top  of  the  bed,  and 
finally  burst  near  the  position  corresponding  to  the  static  bed  height.  The 
small  size  particles  suspended  in  the  air  stream  were  fluctuating  and  weaving 
at  the  top  of  the  bed,  and  so  the  disperse  phenomenon  dominated  the  top  layer 
of  the  bed.  Thus,  there  is  no  clear  height  of  the  fluidized  bed  that  can  be 
measured  visually. 

Three  methods  can.  be  used  to  measure  the  fluidized  bed  height  from  the 
solid  density  profile  under  measurement. 
(a).  End  point  of  the  density  profile: 

Tracing  the  density  profile  in  the  falling  density  zone  and  extrapolating 
to  the  position  at  which  the  density  diminishes  to  zero,  we  can  define  the 
fluidized  bed  height  at  this  point.  This  height  should  correspond  to  the 
maximum  height  of  the  fluidized  bed  as  observed  visually  by  previous  investi- 
gators . 
(b).  Cross  point  of  the  normalized  density  profile: 

As  indicated  in  the  dimensionless  density  profiles,  Figure  13,  the  signi- 
ficant height  of  aggregative  fluidized  beds  can  be  defined  at  the  point  where 
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the  ideally  expanded  density  profile  crosses  the  actually  developed  density 
profile  of  the  fluidized  bed.  The  particles  are  assumed  to  be  uniformly 
distributed  below  this  height,  and  the  portion  above  it  is  less  significant 
to  the  over-all  performance  of  the  total  bed  and  is,  therefore,  considered  to 
be  negligible.  This  height,  as  defined,  should  signify  the  bed,  though  aggre- 
gative, as  a  particulate  one. 
(c).  Height  determined  by  IIS  methods j 

The  fluidized  bed  height  -  defined  as  the  breaking  point  of  the  axial  IIS 
distribution  curve  -  was  proposed  by  Steward  (25).  However  this  height  corres- 
ponds closely  to  the  height  determined  by  method  (b). 

It  is  expected  that  the  bed  defined  by  method  (b)  or  method  (c)  would  give 
better  correlation  than  by  method  (a),  if  the  analogous  correlation  of  a  par- 
ticulate fluidized  bed  is  applied  in  an  aggregative  fluidized  bed.  In  the 
calculation  for  reactor  conversion,  the  bed  height  defined  in  (b)  and  (c) 
should  give  a  more  realistic  picture,  since  the  insignificant  part  of  the 
disperse  phase  is  neglected  in  this  consideration. 

The  fractional  fluidized  bed  expansion  (L  -  Lwp)A  f>  calculated  from  the 
height  as  determined  by  method  (a)  was  plotted  as  a  function  of  superficial  air 
velocity  with  L^fA^c   as  parameters  in  Fig.  26.  This  correlation  shows  a  linear 
relationship,  which  is  consistent  with  the  data  presented  by  Minet  et.al.for 
chars  (17).  The  effect  of  the  static  bed  height  on  the  expansion  of  fluidized 
beds  is  also  shown  in  this  figure.  Similar  correlation  was  attempted  with  the 
bed  height  defined  by  method  (b)  or  method  (c).  The  result  shows  that  it 
deviates  considerably  from  the  linear  relationship  (Fig.  27).  Fig.  28  is  a 
comparison  between  the  fluidized  bed  height  determined  by  method  (a)  and 
method  (c).  It  indicates  that  for  a  certain  value  of  L  -/D  and  particle 
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size,   these  two  defined  heights  are  different  from  each  other  by  a  propor- 
tionality factor. 

Pressure  Gradient  Correlations 

The  pressure-dr op-flow  diagram  in  a  fluidizing  system  is  well  explained  in 

the  following  paragraph  (6): 

If  the  behavior  of  packed  bed  is  followed  when  a  gas  stream 
is   introduced  from  the  bottom,   the  pressure  drop  at  first   increases 
with  flow  rate  without  bed  expansion.     When  the  flow  rate  is 
sufficient  to  cause  a  pressure  gradient  equal  to  the  buoyant 
weight  of  the  solids  per  unit  volume,   further   increase   in  the 
flow  rate  causes  the  bed  to  expand  —   i.e.    increase  the  fractional 
void  volume  —  the  pressure  drop  remains  substantially  constant.... 

The  two  pressure-drop-flow  diagrams   obtained  are  shown  in  Figs.    28-a.    and 
28-b.     They  resemble  in  shape  to  the  pressure-drop-flow  diagram  for   ideally 
fluidizing  systems   (11),    i.e.,  the  pressure  drop  becomes  essentially  just  the 
weight  of  the  bed  at  the  air  velocity  just  beyond  the  minimum  velocity  of 
fluidization.     This  probably  indicates  that  the  fluidizing  bed  developed  with 
the  air  distributor  (canvas  filter  cloth)  was  macroscopically  almost  equivalent 
to  the  ideal  expansion  of  a  packed  bed.   since  the  size  of  air  bubbles  formed 
was  small  and  the  bubbles  were  nicely  distributed  in  the  bed  at  the  air 
velocity  around  the  minimum  velocity  of  fluidization.     This  also   implies  that 
the  tube-bundle  theory  (h)  developed  for  the  study  of  the  pressure  drop 
through  a  packed  bed  may  be  applied  analogously  to  the  analysis   of  the 
pressure  gradient   in  the  present  fluidizing  systems,   provided  the  fluidized 
bed  height  can  be  correctly  determined. 

As  mentioned  previously,  when  the  air  velocity  reaches  the  minimum  velocity 
of  fluidization,   the  pressure  drop,    A  P  is  constant  and  substantially  equal 
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to  the  weight  of  the  particles  constituting  the  bed.  But  the  pressure  gradient 
is  varied  with  increased  air  velocity,  since  the  air  flow  beyond  the  minimum 
fluidizing  velocity  causes  the  bed  to  expand.  Thus,  the  pressure  gradient 
/V  p/t  decreases  with  the  increase  in  the  air  velocity. 

In  Fig.  29,  the  pressure  gradients,  ^P/L,  are  correlated  in  a  form  of 
friction  factor,  f.  ,  originally  proposed  (6)  for  correlating  the  pressure  drops 

•  „m£e   -l 


through  packed  beds,  against  the  modified  Reynolds  number,  N   -  ~ji 
But,  the  numerical  values  of  the  pressure  gradient,  j(\  P/L,  for  fluidized  beds 
should  depend  on  the  definition  of  fluidizing  bed  heights  used.   (See  previous 
section.)  The  correlation  obtained  with  the  height  determined  by  method  (a) 
yields  the  curves  resembling  the  Ergun  equation  for  packed  bed  in  shape  (U), 
and  have  varying  slopes.  Whereas,  the  correlation  with  that  determined  by 
method  (b),  like  Blake-Kozeney  equation,  has  a  constant  slope.  In  both 
cases  the  observation  of  Fig.  29  and  the  data  shows  that  the  effect  of  particle 
size,  D  ,  appears  to  be  more  significant  than  that  suggested  by  the  correlations 
proposed  for  packed  beds.  The  experimental  data  indicate  that  the  static  bed 
height  also  has  some  effect  on  the  correlation  of  pressure  gradients.  But, 
this  factor  does  not  come  into  the  consideration  for  packed  beds. 

The  results  of  the  dimensionless  analysis  lead  to  the  following  form 
(Appendix  II): 

%  .*£p .    O  .  ^  -  ^  &££  .  _i_  )a&)b(^)c 

*    L     1-6         PpV2  T      S*  l-£     Dc   Dc      (2h) 

Where,  the  constant  (p   ,   a,  b,  c,  can  be  evaluated  from  the  experimental  data 
with  the  height  of  fluidized  bed  defined  by  methods  (a)  and  (b). 

With  the  height  defined  by  method  (a),  the  dimensional  analysis  gives: 
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k  -  —  •  rrr '  7^ ~ 

S    Dd/Vv  1         -°'^    °t>    1.389^  -0.22U 

(9.028  x  10*)(^L  .  r^-)  (£)  (-«£)  (25) 

and,  with  the  height  defined  by  method  (b): 

9rAP         ^  ^  D« 

K         L  l-£        jOpV2 

(8.80  x  H>')(-£^-  .  -i_)  &L-)  (JS£)  (26) 

Equations  (25)  and  (26)  are  plotted  in  Figs.  (30)  and  (31)  respectively. 

The  comparison  plot  for  Equations  (25)  and  26)  is  shown  in  Fig.  32  in 
order  to  see  which  definition  of  height  for  fluidized  bed  gives  better  results 
in  the  proposed  correlation  for  the  pressure  gradient.  It  can  be  seen  from 
this  plot  that  the  height  determined  by  method  (a)  gives  more  scattered 
correlation  than  that  determined  by  method  (b).  It  may  be  concluded,  in  the 
proposed  dimensionless  correlation  of  pressure  gradient  through  fluidized 
bed,  that  the  height  determined  by  method  (b)  gives  more  consistent  result 
than  that  by  method  (a). 
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RESULTS  AMD  DISCUSSION 


£T-ray  Technique  as  Applied  to  the  Measurement 
of  Fluidized  Bulk  Density  in  Fluidized  Beds 


Among  the  several  technical  methods  developed  for  the  study  of  local 
properties  of  a  non-ideal  two  phase  fluid,   the     ^-ray  technique  is  a  rela- 
tively new  one.     It  could  be  applied  to  the  study  of  the  formation  and 
behavior  of  bubbles  in  aggregative  fluidized  beds.     By  using  a  short-time 
constant  of  the  rate-count-meter  to  obtain  the  strip  chart  oscillograms, 
the  size,   velocity  and  frequency  of  bubbles  could  be  calculated  (3).     This 
method  also  simplified  the  technical  problems   in  the  measurement  of  local 
time-average  density  by  using  a  long-time  constant.     The  capacitance  probe 
method  which  was  used  by  Dotson  (5)  and  other  investigators   (18,   1.) 
appeared  to  disturb  the  flow  pattern  of  the  fluidizing  vessel.     The  capaci- 
tance probe  inserted  in  the  vessel,   though  it  was  small,   should  change  the 
local  flow  pattern.     The  properties   investigated  were  influenced  by  the 
probe  and  should  be  different  from  the  flow  pattern  without  such  disturbance. 
The    "25"-ray  radiation  technique  can  get  rid  of  this  disturbance  completely. 
The  radiation  passing  through  the  fluidizing  vessel  has  no  effect  on  the 
motion  of  fludizing  mixtures. 


Comparison  of  the  Short-time  Constant 
to  Long-time  Constant  Data 


Although  a  30  second  long  time  constant  was  used  for  all  the  experimental 
runs,   the  data  obtained  were  compared  with  the  data  taken  by  short-time 
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constant   (0.3  second)   in  one  particular  run.     The  method  used  to  calculate 
the  line  average  density  from  the  fluctuating  CRM  reading  with  a  short-time 
constant  was  essentially  the  same  as  that  used  by  Stewart   (25).     The  axial 
density  profiles  obtained  with  both  long-time  constants  and  short-time 
constants  are  plotted  in  Fig.   3h  for  comparison.     It  can  be  seen  that  the 
results  were  almost  identical,    indicating  that  the  sensitivity  of  the  rate- 
count-meter  was  satisfactory.     Considering  the  effort  required  to  determine 
the  axial  density  profile  from  the  short-time  constant  data,   we  should  be 
satisfied  with  the  results   obtained  with  a  long-time  constant  for  the 
purpose  of  the  present  studies.     The  density  profiles  were  most  directly 
obtained  from  the  recording  chart   in  long-time  constant  runs. 

Reproducibility 

The  accuracy  of  the  experimental  result  was  determined  by  the  consis- 
tency test.     The  consistency  test  was  satisfactory  if  the   instrument  operated 
was  stable  and  the  power  to   detect  the   ^f-ray  was  consistent  throughout. 
Occasionally,   the  instrument  became  damped  during  the  process  of  experiment, 
and  the  density  profile  obtained  would  be  quite  skewed  even  in  the  constant 
density  zone.     This  is  shown  in  Fig.   35*     This  deficiency  could  be   improved 
by  checking  the  calibration  before  and  after  taking  the  data,    and  by  avoiding 
any  mechanical  shock  to  the    y~ray  source.     In  fact,   the  reproducibility 
obtained  was  good,   considering  that  the  fluidizing  state  of  beds  changed 
constantly  owing  to  the  nature  of  bed  itself.     Installing  more  than  one 
2T~ray  source  and  scintillation  probe  to  measure  the  solid  density  from 
several  directions  of  the  fluidizing  column  would  have  been  necessary  to   im- 
prove the  reproducibility  to  any  larger  extent. 
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The  Importance  of  the  "Bed  Density  Distribution" 

to  the  Understanding  of  Heat  and  Mass  Transfer  Mechanism 

in  Catalyst  Fluidized  Bed  Reactor 


In  spite  of  the  wide  application  and  rapid  development  of  fluidization 
techniques,  a  complete  understanding  of  the  fundamental  factors  governing  the 
performance  of  fluidizing  units  has  been  lacking  (especially  for  gas-solid 
systems.)*  One  way  for  expressing  the  bed  performance  is  in  terms  of  various 
degrees  of  so-called  "good"  or  "bad"  fluidization  based  on  some  definite  index- 
This  index  is  a  direct  measure  of  the  effectiveness  of  heat  and  mass  transfer 
as  well  as  the  conversion  of  reactants  in  chemical  reactors  (l£,.l6,  18), 

Heat  and  mass  transfer  correlations  for  packed  beds  are  rather  complete 
in  the  literatures  (11.,  28),  There  has  been  a  trend  during  the  past  15  years 
for  investigators  in  the  field  of  fluidization  to  use  the  analogy  to  the  packea 
bed  reactor  in  the  study  of  fluidized  bed  reactors.  But  the  difficulties 
encountered  were  due  to  the  non-uniform  particle  distribution  throughout  the 
fluidized  bed.  The  reacting  gas  flowing  through  the  fluidized  bed  is  not  only 
a  tortuous  flow  through  the  constant  void  space  as  in  packed  beds,  but  in  addi- 
tion it  is  affected  by  the  impacts  between  particles  and  the  interaction  between 
particles  and  the  flowing  medium.  A  quantitative  measurement  of  the  complexity 
of  the  flowing  pattern  in  the  fluidized  bed  is  therefore  important.  The  mass 
and  heat  transfer  in  the  fluidized  bed  is  governed  by  the  surface  and  the  move- 
ment of  particles  in  the  bed.  It  is  certain  that  the  transfer  is  better  in 
the  locations  where  the  solid  particles  are  dense  and  are  constantly  circu- 
lating. 

Mickley  et.  al.  (16)  in  their  study  of  heat  transfer  for  gas  fluidized  beds 
concludes: 
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The  local  time-mean  wall-to-bed  heat  transfer  coefficient  varies 
with  the  axial  position  of  heat  transfer  surface,   decreasing  with 
increasing  distance  above  the  bottom  of  the  bed.     On  the  other  hand, 
the  local  coefficient  is  not  significantly  influenced  by  total  height 
of  the  bed. 

Shirai  et»al.(23)  also  drew  the  same  conclusion  in  studying  the  axial 
distribution  of  heat  transfer  coefficients.     The  axial  distribution  of 
particles-to-fluid  heat  transfer  coefficients   obtained  by  them  is  almost 
identical  in  shape  to  the  axial  density  profile  obtained  in  this  work. 

The  close  relationship  between  the  axial  distribution  of  the  heat 
transfer  coefficient  and  the  axial  density  profile  suggests  that  the  study 
of  the  density  profile  of  fluidized  beets  will  help  in  understanding  heat 
and  mass  transfer  mechanism  in  the  fluidized  bed,     since  the  particle  behavior 
is  the  major  factor  affecting  the  performance  of  fluidized  beds. 


CONCLUSION 

The  following  significant  conclusions  have  been  drawn  from  the  results 
of  the  present  investigations 

1.  With  a  long-time  constant  on  the  CRM,  the     ^T-ray  technique  can  be 
applied  satisfactorily  to  the  study  of  the  axial  density  profile  of  the 
aggregative  fluidized  bed.     The  reproducibility  is  within  the  allowable  error, 
except  those  at  very  high  air  velocities. 

2.  The  data  taken  with  a  short-time  constant   is  almost  identical  with 
the  data  using  long-time  constant. 

3.  Two  distinct  zones  can  be  seen  in  the  axial  density  profile,   the 
constant  density  zone  and  the  falling  density  zone.     The  line  average  density 
in  the  constant  density  zone  decreases  with  increasing  air  velocity. 


h>     The  ratio  of  the  static  bed  height  to  the  column  diameter   instead  of 
the  static  bed  height   itself   is   the  important  variable   influencing  the 
fluidizing  state. 

5.  Indices  A  and  B  can  be  used  to  describe  aggregative  fluidization. 
The  greater  the  Index  B,    the  more  stable   is  the  fluidizing  bed. 

6.  Indices  A  and  B  are  correlated  by  the  following  empirical  equations: 

I  A  I  -  1.3U9  x  10-5(«     )0.978(i)-l.a6S(W)0.3W9  (22) 

^C  ^c 

I  B  I  -  5.U92  x  io9(NRe)-2-0^(^)3-OW(W)0.3l17l»  {2}) 

7.  Among  the  three  significant  operational  variables,  air  velocity, 
static  bed  height,  and  particle  size,  air  velocity  has  the  most  pronounced 
effect  on  the  fluidizing  state. 

8.  The  bed  height  defined  as  the  end  point  of  the  axial  density  profile 
corresponds  to  the  visual  maximum  bed  height.  The  more  significant  mean  height 
of  aggregative  fluidized  bed  is  defined  as  the  cross  point  of  the  normalized, 
density  profile. 

9.  Pressure  gradient  (^P/L)  correlation  yields  the  following  equations? 

f„-3>      9£££._£_L- 

K       j>2         L  l-£ 

S     DdvPf  1         -0.9W  Dn  1.389  .-  -0.22li 

(9.028  x  105)(-^P  '  j^)  (-E)  q*)  (25) 

m    Dp       1 9cAP  _     £3 

K  =  fFv*     l    '  1  -e 

8     DpVpF  !       -1.591  dd  2.126  Lmf  -0.298 

(8.80  x  lo8)(-^p  '  j-Jg-)  (^)  (JE)  (26) 

In  equation  (25),  L  was  determined  by  method  (a),  and  in  equation  (26),  by 
method  (b)  or  (c).  The  second  equation  correlates  the  data  better  than  the 
first  one. 


91 


10.     The  heterogeneous  distribution  of  fluidized  bulk  density  in  the 
aggregative  fluidized  beds  has  a  direct  influence  on  the  heat  transfer,  mass 
transfer  and  local  conversion  of  fluidized  bed  reactors. 
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NOMENCLATURE 

o 

A  =  cross-sectional  area  of  the  fluidized  bed,  1  . 

p 

A'  =  part  of  cross-sectional  area  of  the  fluidized  bed,  1  . 

Ar  «  projected  area  of  ^T-beara  path,  1  . 

A0  =  constant,  dimensionless. 

B0  ■  constant,  dimensionless. 

A  =  index  of  fluidized  bed,  dimensionless. 

B  ■  index  of  fluidized  bed,  dimensionless. 

d^^  s  thickness  of  spacer  aluminum,  1. 

Dc  =  inner  diameter  of  fluidized  column,  1. 

dc  =  wall  thickness  of  fluidized  column,  1. 

D_  =  diameter  of  particle,  1. 

Jr 

Dp'  ■  width  of  packed  particle  bed  in  the  direction  of  radiation  beam,  1, 

Ds  ■  width  of  spacer,  1. 

fjr  =  friction  factor,  dimensionless. 

G  «  volumetric  air  velocity,  1-7©  • 

G,  s  volumetric  bubble  velocity,  l-y9  . 

g  -  conversion  from  force  to  mass,  ml/F0  . 

G^  B  volumetric  velocity  of  minimum  fluidization,  1^/0  . 

h  «  height  above  the  distribution  in  bed,  1* 

l/l0  ■»'  fraction  of  photons  remaining  in  the  beam  after  passage  through 
absorber  of  thickness  d,   dimensionless. 

L  «  height  of  fluidized  bed,  1. 

L  -  ■  static  (on-set)  height  of  bed,  1. 

M  =  total  particle  weight  in  bed,  m. 

M'  =  bed  weight  obtained  from  integrated  area  (consistency  test),  m. 
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Ms  a  weight  of  particles  fluidizing  above  L^f,  m. 

KL  «  weight  of  particles  fluidizing  below  Lp_,  m. 

Na  «  Avagodro*s  number   (6.025  x  10  ^  nuclei/g-mol). 

NR  e  Reynolds  number,  V  PpD_/^/      ,   dimensionless. 

%e  ""  m°dified  Reynolds  number,  vPpD  A^        *   l/l  -£        ,   dimensionless. 

Z^   P  ■  pressure  drop  through  fluidized  bed,   F/l2 

t^  «  the  residence  time  of  bubbles   in  the  bed  below  L  f,   9   . 

V  s  superficial  air  velocity,   1/0. 

V]-,  ■  average  linear  velocity  of  bubbles,   I/O  • 

Vinf  a  minimum  velocity  of  fluidization,    l/O. 

Vol]-,  ■  the  volume  occupied  by  bubbles  below  L^,    1   . 

Vol  ■  volume  of  particles  fluidizing  above  L^,   1^. 

Vols  a  volume  of  particles  fluidizing  below  L^f,    1   • 


Greek  Letters 


C^ »  -  constant . 

p  t  a  constant. 

£    ~  a  porosity  of  the  fluidized  bed  at  on-set  fluidization. 

P  B  fluidized  bulk  density. 

P,»  ■  density  of  aluminum. 

P  «  density  of  the  fluidizing  column  material. 

P  »  fluidized  bulk  density  at  the  position  corresponding  to  the  static 

'  bed  height. 

Pd  "  average  line  density  in  the  constant  density  zone. 

Pp  «  density  of  fluid. 

r^.  ■  fluidized  bulk  density  at  on-set  fluidization. 

P   a  «  bulk  density  of  bed  at  packed  state. 
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P  «  density  of  solid  particle. 

(D  »  constant. 

u  h  microscopic  absorption  cross-section. 

S*  »  viscosity  of  fluid. 

>^/  m  ~  absorption  coefficient. 
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APPENDIX  II 

(A)     Derivation  of  the  Dimensionless  Equation  for  Correlating  Pressure 
Gradient   in  the  Aggregative  Fluidized  Bed -Friction  Factor,   fK 

The  pressure  gradient,    Z^  P/L,   has  been  found  experimentally  to  be  in- 
fluenced by  the  following  properties   of  column,   fluid,   and  particles s     column 
diameter,  Dc,   fluid  density,   Pp,    viscosity,  Li  ,   velocity,  V,  particle  diameter, 
D„,   and  static  bed  height,   L^f .     The  primary  relationship  existing  between  the 
pressure  gradient  and  the  variables   is   (6,    11,    28) % 

Ap        ,  h    n     {  jkm  np 

-^  «  <£(DC)   (/F)    (JU.)   00    (Dp)    (Lmf)    (gc)  (27) 

where    (P  ,   h,    i,    j,   k,   m,   n,   q,   are  unknown  constants,   and  gc   is  the  con- 
version constant  between  force  and  mass. 

Equation  (27)   is  rewritten  using  the  dimensions  of  each  variable  ass 

hi  j  <  k      m      n     .,     q 

^  -4>  (i)   (f)   (y-jj^)    P)   (1)   (1)    (~s)  (28) 

■Summing  the  exponents  of  equation  (28), 

IF,  1  -  -  q 

2  1,  »3~h-3i-j  +  k  +  m  +  n  +  q 

2Tm,  0  =  i  +  j  +  q 

2  0,  0--J-k-2q 

To  solve  the  above  four   equations  for  seven  unknowns,   the  Bridgman 
Pi-Theorem'"  is  used.     Choosing  h,   j,   n,   as  three  unknowns,   then, 


'Bridgman's  Pi-theoremj  If  there  are  N  number  of  equations  to  solve  for  K 
unknowns,  the  number  of  dimensionless  groups  to  be  obtained  is  W  -  K  +  1, 
and  M  -  K  of  unknowns  retain     numerically  unsolved. 


q  =  -   1 

iB-  j-  qB  1-  j 
k--j-2q-2-J 
m°-3-h  +  3i  +  j-k-n-q 
a-h-j-n-1 
Substituting  the  above  values   into  equation  (27)s 


ap  .         h  1-j        J       2-j        -1-h-j-n  n        -1 

(^L£)  «  0(Dc)   (j>p)       (Ji)  (V)       (Dp)  (Lmf)  (gc) 


and  rearranging: 


or, 


f3cAP  ,  "p     a  f-PFVDp.-JPp.-h-n.Lmf.n 


Let  -  j  =  a,  -  h  -  n  *  b,  n  =  c,  then: 

a 


j>Fv2 '     r  ^ 


c   ^c 


10 


? 


(29) 


go^_P  .  _Dp_        ,  -ZaSe,  A)  (hi)0  00) 


Based  on  the  theoretical  consideration  (for  the  hydraulic  radius  of 
column  and  the  wetted  surface  of  particles,  etc.)>  the  friction  factor,  f^, 
derived  for  the  packed  bed  has  been  modified  by  multiplying  the  dimensionless 
factor,  £  V1  -  £   (6),  and  the  Reynolds  number,  'PpVDp/^.  ,  has  been 
modified  by  multiplying  by  l/l  -  6  (where  £  is  the  void  fraction  of  the 
fluid-particles  bed.)  The  dimensionless  equation  then  becomes: 
9c^P   Dp     £3 

-  ^(J^2.   1   )m(3»b££)0  (ft) 

7   •"     1-6    Dc   Dr 


'c 


log 


(B)  Derivation  of  the  Dimension!. ess  equation  for  Correlating  Indices  A  and 
B. 

The  results  of  the  previous  and  the  present  investigation  indicate  that 

(5,  25): 

h     i     j   k    m     n 
|  A|  -  0(  «(DC)  (j°F)  {/L  )  (V)  (Dp)  (Lmf)  (31) 

h    5  i        j     k   m   n 
/  -  Ol  '(0  (VlJ)  (»/l  -  0)  (VO)  (1)  (1) 

5"  m,  0  ■  i  +  j 

2"  1,  0  =  h-3i  -  j  +  k  +  m  +  n 

10,   0  -  -  j ■  -  k 
Holding  h.,    j,    n,    as  unkno^m  constants 

i  =  -  j 

k  -  -  j 

m  -  -  h  +  3±  +  j  -  k  «  n  =  -  h  -  j  -  n 
substituting  these  into  equation  (31), 

|A|   -  d   ,(Dcy1(/F)"J(/^)J(V)"J(DF) J  "  n(Lmf)n  (32) 

Rearranging: 

I  A  |       -  o<   » (//  If  FVDp)  J (Dc /Dp)    (Lmf /Dp )n 

«(*    «(/FVDp/^      )"   (Dp/Dc)'  n(Lmf/Dc)n 

Let  -  j  ■  a,   -  h  -  n  =  b,  n  -  c,   then: 

I A  |      -cX  '(/FVDp/^      )    (Dp/Dc)    (Lmf/Dcf  (20) 

By  the  similar  method,   the  dimensionless  equation  for  correlating   index  B  can 
be  obtained  as: 

IB  |      "    £   '(jVV>     )a'(VDc)b'(Lmf/Dc)C'  <21> 
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q   -ray  attenuation  method  was  applied  to  determine  the  axial  fluidizing 
bulk  density  profile  for  an  air-solid  aggregative  fluidized  bed.     Two  dis- 
tinct density  zones,   constant  and  falling,  were  observed  except  for  the  bed 
fluidizing  at  very  high  air  velocity. 

The  average  density  of  the  constant  zone    P^,  was  correlated  as  a  func- 
tion of  the  operating  variables  ass 

f  d  a     01-/3     (V  "  Vmf ) 

where  0(  and  u  are  the  constants  depending  on  the  particle  diameter  only, 
V  and  Vflf  are  fluidizing  gas  velocity  and  minimum  fluidizing  velocity  of  the 
particles. 

The  density,  P  ,  in  the  falling  zone  was  found  to  be  approximated  by 
the  following  equation,  -which  is  characterized  by  the  two  indices  A  and  Bs 

2 


Log  (f/f   pa)  =  A  +  B  (h/L^  -  1)' 


The  effects  of  operational  variables,  such  as  air  velocity,  V,  static  bed 
height,  L~,  particle  size,  D  ,  on  these  indices  were  investigated.  They 
were  expressed  as  functions  of  dimensionless  groups  ass 
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The  effects  of  particle  composition  in  beds  containing  a  mixture  of  two 
sizes  of  particles   (1|0~U5  mesh  and  80-100  mesh  glass  beads)   on  fluidization 
characteristics  were  also  studied.     The  results  showed  that  the  values  of 


indices  A  and  B  were  significantly  changed  only  when  the  mixture  contained 
less  than  10  weight  percent  of  one  component. 

Three  different  ways  to  define  the  height  of  a  fluidized  bed  were  con- 
sider eds 

(a)  End  point  of  the  density  profile. 

(b)  Cross  point  of  the  normalized  density  profile  with  that  of  an 
ideal  fluidized  bed. 

(c)  Maximum  point  of  axial  instability  index  distribution  curve. 

The  height  defined  in  (a)  corresponded  to  the  maximum  bed  height  as  visually 
observed,  while  the  other  two  are  the  average  of  the  fluctuating  bed  height. 
Pressure  gradient  was  correlated  in  the  form  of  a  friction  factor.  The 
heights  defined  in  (b)  and  (c)  gave  better  correlation  than  that  defined  in 
(a).  The  friction  factor  obtained  is  given  as  follows  % 
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The  heterogeneous  distribution  of  fluidized  bulk  density  in  the  aggre- 
gative fluidized  beds  has  direct   influences   on  the  transport  and  rate 
processes.     Therefore,  the  knowledge  of  the  density  distribution  in  the 
aggregative  fluidized  beds  cannot  be  over -emphasized.     The  works  of  previous 
investigators  on  the  transport  and  rate  processes   in  fluidized  beds  were 
re-examined  in  the  light  of  the  results  obtained  in  this   investigation. 


